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Lower Bounds on Complexity of Lyapunov 
Functions for Switched Linear Systems 
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Abstract 

We show that for any positive integer d, there are families of switched linear systems— 
in fixed dimension and defined by two matrices only—that are stable nnder arbitrary 
switching but do not admit (i) a polynomial Lyapunov function of degree < d, or (ii) a 
polytopic Lyapunov function with < d facets, or (iii) a piecewise quadratic Lyapunov 
function with < d pieces. This implies that there cannot be an upper bound on the 
size of the linear and semidefinite programs that search for such stability certificates. 
Several constructive and non-constructive arguments are presented which connect our 
problem to known (and rather classical) results in the literature regarding the finiteness 
conjecture, undecidability, and non-algebraicity of the joint spectral radius. 

In particular, we show that existence of an extremal piecewise algebraic Lyapunov 
function implies the finiteness property of the optimal product, generalizing a result 
of Lagarias and Wang. As a corollary, we prove that the finiteness property holds for 
sets of matrices with an extremal Lyapunov function belonging to some of the most 
popular function classes in controls. 

Index terms: stability of switched systems, linear difference inclusions, the finiteness 
conjecture of the joint spectral radius, convex optimization for Lyapunov analysis. 


1 Introduction 

In recent years, hybrid systems have been well recognized as a powerful modeling 
framework for complex dynamical systems. This in turn has led to a surge of 
activity to design computational tools for their analysis and control. Many 
of these tools are developed based on a synergy between classical ideas from 
Lyapunov theory (often appropriately modified or generalized to deal with the 
intricacies of hybrid systems) and modern techniques of numerical optimization, 
notably convex optimization. The outcome of this synergy is an automated 
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search process for appropriately-defined Lyapunov functions that certify desired 
properties of system trajectories (e.g., stability, boundedness, etc.). 

Our goal in this paper is to show some basic and concrete limitations of this 
broad class of techniques for analysis of hybrid systems. Since we are aiming at 
negative results, the simpler we can make the mathematical setting of our study, 
the more powerful the implications of our theorems will be. For this reason, 
we choose to focus on arguably the simplest class of hybrid systems, namely, 
switched linear systems, and ask arguably the simplest question of interest about 
them, namely, stability. 

More precisely, the input to our problem is a set of m real n x n matrices 
S := {Ai ,..., Am}. This set describes a switched linear system of the form 

— A(j(^]^'jX}^, ( 1 ) 

where k is the index of time and a : Z —>-{1, ...,m} is a map from the set of 
positive integers to the set {!,...,m}. A basic notion of stability is that of 
absolutely asymptotically stable (AAS), also referred to asymptotic stability un¬ 
der arbitrary switching (ASUAS), which asks whether all initial conditions in 
R" converge to the origin for all possible switching sequences. It is not diffi¬ 
cult to show that absolute asymptotic stability of © is equivalent to absolute 
asymptotic stability of the linear difference inclusion 

Xk+i G coS Xk, (2) 

where coS here denotes the convex hull of the set S. Dynamical systems of 
this type have been widely studied since they model a linear system which is 
subject to time-dependent uncertainty. See for instance [37] j EHIj or [3T] for 
applications in systems and control. 

When the set S consists of a single matrix A (i.e., m = 1), we are of course 
in the simple case of a linear system where asymptotic stability is equivalent to 
the spectral radius of A having modulus less than one. This condition is also 
equivalent to existence of a quadratic Lyapunov function and can be checked 
in polynomial time. When m > 2, however, no efficiently checkable criterion 
is known for AAS. Arguably, the most promising approaches in the literature 
have been to use convex optimization (typically linear programming (LP) or 
semidefinite programming (SDP)) to construct Lyapunov functions that serve 
as certificates of stability. The most basic example is that of a common quadratic 
Lyapunov function (CQLF), which is a positive definite quadratic form x"^Qx 
that decreases with respect to all m matrices, i.e., satisfies 

x^{AfQAi — Q)x < 0,Va; G ]R"'\{0},i = 1,... ,m. 

On the positive side, the search for such a quadratic function is efficient numeri¬ 
cally as it readily provides a semidefinite program. On the negative side, and in 
contrast to the case of linear systems, existence of a CQLF is a sufficient but not 
necessary condition for stability. Indeed, a number of authors have constructed 
examples of AAS switched systems which do not admit a CQLF and studied 
various criteria for existence of a CQLF ( [71II11IM1I3I])- 
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To remedy this shortcoming, several richer and more complex classes of 
Lyapunov functions have been introduced. We list here the five that are perhaps 
the most ubiquitous: 

Polynomial Lyapunov functions. A homogeneous0 multivariate polyno¬ 
mial p{x) of some even degree d is a polynomial Lyapunov function for ([1]) if it 
is positive definit^l and makes p{x) — p{Aix) positive definite for t = 1,..., m. 

Although this is a rich class of functions, a numerical search for polynomial 
Lyapunov functions is an intractable task even when the degree d is fixed to 
4. In fact, even testing if a given quartic form is positive definite is NP-hard 
in the strong sense (see, e.g., [T]). A popular and more tractable subclass of 
polynomial Lyapunov functions is that of sum of squares Lyapunov functions. 

Sum of squares (sos) Lyapunov functions. A homogeneous polynomial 
of some even degree d is a sum of squares (sos) Lyapunov function for o 
if p is positive definite and a sum of square 4 j 7 and if all the m polynomials 
p{x) — p(Aix), i = 1,... ,m are also positive definite and sums of squares. 

For any fixed degree d, the search for an sos Lyapunov function of degree 
d is a semidefinite program of size polynomial in the input. When d = 2, this 
class coincides with CQLFs as nonnegative quadratic forms are always sums of 
squares. Moreover, existence of an sos Lyapunov function is not only sufficient 
but also necessary for AAS of o m- This of course implies that existence 
of polynomial Lyapunov functions defined above is also necessary and sufficient 
for stability. 

Polytopic Lyapunov functions. A polytopic Lyapunov function V for 
CD with d pieces is one that is a pointwise maximum of d linear functions: 

V{x) := max \cfxl 

where ci,..., Cm span R”. The sublevel sets of such functions are polytopes, 
justifying their name. Polytopic Lyapunov functions (with enough number of 
pieces) are also necessary and sufficient for absolute asymptotic stability. One 
can use linear programming to search for subclasses of these Lyapunov functions. 
These subclasses are big enough to also comprise a necessary and sufficient 
condition for stability (see [51 I^I?41I55] 1. 

Max-of-quadratics Lyapunov functions. A max-of-quadratics Lyapunov 
function V for CD with d pieces is one that is a pointwise maximum of d positive 
definite quadratics: 

V{x) := max x'^QiX, 

where Qi >- 0. The sublevel sets of such functions are intersections of ellipsoids. 

^ Since the dynamics in 0 is homogeneous, there is no loss of generality in parameterizing 
our Lyapunov functions as homogeneous functions. Also, we drop the prefix “common” from 
the terminology “common polynomial Lyapunov function” as it is implicit that our Lyapunov 
functions are always common to all m matrices Ai in S. 

^ A form (i.e., homogeneous polynomial) p is positive definite if p{x) > 0 for all x ^ 0. 

^ A polynomial p is a sum of squares if it can be written as p = some polynomials 
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Min-of-quadratics Lyapunov functions. A min-of-quadratics Lyapunov 
function V for o with d pieces is one that is a pointwise minimum of d positive 
definite quadratics: 

V(x) := min x'^QiX, 

where Qi >- 0. The sublevel sets of such functions are unions of ellipsoids. 

By a piecewise quadratic Lyapunov function, we mean one that is either 
a max-of-quadratics or a min-of-quadratics. Both of these families are known 
to provide necessary and sufficient conditions for AAS. Several references in 
the literature produce semidehnite programs that can search over a subclass of 
max-of-quadratics or min-of-quadratics Lyapunov functions (see [17jl. These 
subclasses alone also provide necessary and sufficient conditions for AAS. A 
unified framework to produce such SDPs is presented in 0, i, where a recipe 
for writing down stability proving linear matrix inequalities is presented based 
on some connections to automata theory. 

For all classes of functions we presented, one can think of d as a complexity 
parameter of the Lyapunov functions. The larger the parameter d, the more 
complex our Lyapunov function would look like and the bigger the size of an 
LP or an SDP searching for it would need to be. 

1.1 Motivation and contributions 

Despite the encouraging fact that all five classes of Lyapunov functions men¬ 
tioned above provide necessary and sufficient conditions for AAS of H]) that are 
amenable to computational search via LP or SDP, all methods offer an infinite 
hierarchy of algorithms, for increasing values of d, leaving unclear the natural 
questions: How high should one go in the hierarchy to obtain a proof of stabil¬ 
ity? How does this number depend on n (the dimension) and m (the number 
of matrices)? Unlike the case of CQLF which is ruled out as a necessary condi¬ 
tion for stability through several counterexamples in the literature, we are not 
aware of that many counterexamples that rule out more complicated Lyapunov 
functions. For example, is there an example of a set of matrices that is AAS 
but does not admit a polynomial Lyapunov function of degree 4, or 6, or 2OO‘0 
Or, is there an example of a set of matrices that is AAS but does not admit a 
piecewise quadratic Lyapunov function with 200 pieces? If such sets of matrices 
exist, how complicated do they look like? How many matrices should they have 
and in what dimensions should they appear? 

In this paper we give an answer to these questions, providing constructive 
and non-constructive arguments for existence of “families of very bad matrices”, 
i.e., those forcing the complexity parameter d of all Lyapunov functions to be 
arbitrarily large, even for fixed n and m (in fact, even for the minimal nontrivial 
situation n = m = 2). The formal statement is given in Theorem [T] below. 

The largest degree existing counterexample that we know of is one of our own, appearing 
in [2], which is a pair of AAS 2x2 matrices with no polynomial Lyapunov function of degree 
14. 
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It is important to remark that the families of matrices we present have 
already appeared in rather well-established literature, though for different pur¬ 
poses. These matrices have to do with the “non-algebraicity” and the “finiteness 
property” of the notion oi joint spectral radius (JSR) (see Sections H] and [3] for 
definitions). This leaves us with the much simplified task of establishing a for¬ 
mal connection between these two concepts and that of complexity of Lyapunov 
functions. We hope that clarifying these connections sheds new light on the 
intrinsic relationship between the JSR and the stability question for switched 
linear systems. Indeed, many of the results that we refer to in the literature 
on the JSR appear much before counterexamples to existence of CQLF in the 
switched systems literature. 

The theorem stated below is proven over the next two sections. 


Theorem 1. For any positive integer d, the following families of matrices (pa¬ 
rameterized by k) include switched systems that are asymptotically stable under 
arbitrary switching but do not admit (i) a polynomial (hence sos) Lyapunov 
function of degree < d, or (ii) a polytopic Lyapunov function with < d facets, 
or (Hi) a max-of-quadratics Lyapunov function with < d pieces, or (iv) a min- 
of-quadratics Lyapunov function with < d pieces: 

1. (1 - l){Ai,A 2 }, with 


(l-t*) -t 

(I-37rt3/2)[ 0 0_ 


^2 = (1 - t *) 


VT^ 

t 


— t 


where t = sin o,nd fc = 1, 2,.... 

(This family appears in the work of \25^ as an example demonstrating that 
the joint spectral radius is not a semialgebraic quantity; see Section [H ) 


2 . 


with 


Ai 


1 

0 


,A2 — I 


0 

1 ’ 


a < a* ~ 0.749326546330367557943961948091344672091.... 


(This set appears in U^\20f as a counterexample to the finiteness prop¬ 
erty.) 


3. (1 - j){Ai,A 2 }, with 


Ai=a^ 


0 0 
1 0 


, ^2 = a 


-1 


1 < a < (cos—) ^ 


where 
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(This family appears in the work of 126}/ as an example demonstrating that 
the length of the optimal product cannot be bounded; see Section\^) 

The first construction and its relation to non-algebraicity is presented in Sec¬ 
tion [21 The second and third constructions are similar and their relations to the 
finiteness property are presented in Section [31 In Section [31 we present a result 
that is of potential interest independent of the above theorem: that existence of 
an extremal piecewise algebraic Lyapunov function implies the finiteness prop¬ 
erty of optimal products. Piecewise algebraic Lyapunov functions constitute a 
broad class of functions that includes all classes stated so far as a special case. 
By an extremal Lyapunov function, we mean a Lyapunov function V such that 

V{A,,x) 

sup = P, 

where p is the JSR. See mi Section 2.1] for more on the similar concept of 
extremal norms. 

In a sense, this result links lower and upper bound approaches for computa¬ 
tion of the joint spectral radius. Similar results were obtained in the pioneering 
works of m for polytopic Lyapunov functions and [2S] for quadratic Lyapunov 
functions, as well as, several other classes of convex Lyapunov functions. 

We shall also remark that for continuous time switched linear systems, m 
have established that the degree of a polynomial Lyapunov function for an ASS 
system may be arbitrarily high, answering a question raised by Dayawansa and 
Martin. We have been unable to come up with a transformation from continuous 
time to discrete time that preserves both AAS and non-existence of polynomial 
Lyapunov functions of any desired degree. 

In Section [31 we provide an alternative proof of Theorem [1] based on an un¬ 
decidability results due to m- While this will be a non-constructive argument, 
its implications will be stronger. Indeed, Theorem |T] above implies that the 
complexity parameter d (and hence the size of the underlying LPs and SDPs) 
cannot be upper bounded as a function of n and m only. The undecidability 
results, however, imply that d cannot be upper bounded even as a (computable) 
function of n, m, and the entries of the input matrices. 

In Section jS] we numerically solve semidefinite programs to demonstrate the 
result of Theorem [T] with some concrete numbers. It is shown that as some 
parameter in the input matrices changes, the minimum required degree of a 
polynomial Lyapunov function indeed grows rapidly (all the way to 24 before 
we stop the experiment). We close our paper with some brief concluding remarks 
in Section [31 

2 Complexity of Lyapunov functions and non-algebraicity 

One classical approach to demonstrate that a problem is hard is to establish 
that there is no algebraic criterion for testing the property under consideration. 
This is formalized by showing that the set of instances of a given size that satisfy 
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the property do not form a semialgebraic set (see formal definition below). Such 
a result rules out the possibility of any characterization of the property at hand 
that only involves operations on the input data that include combinations of 
arithmetical operations (additions, subtractions, multiplications, and divisions), 
logical operations (“and” and “or”), and sign test operations (equal to, greater 
than, greater than or equal to,...); see m- While this is a very strong statement, 
non-algebraicity does not imply (but is implied by) Turing undecidability, which 
will be our focus in Section S) Nevertheless, non-algebraicity results alone are 
enough to show that the complexity of commonly used Lyapunov functions 
for switched linear systems cannot be bounded. The goal of this section is to 
formalize this argument. 

Definition 1. A set S C K” defined as S = {x : ffix) 0, f = 1,..., r}, 
where for each i, fi is a polynomial and l>i is one of>, <, =, is called a basic 
semialgebraic set. A set is called semialgebraic if it can be expressed as a finite 
union of basic semialgebraic sets. 

Theorem 2 ( [41], [38]). Let S C be a semialgebraic set and tt : — >■ 

R." be a projection map that sends {x,y) i—f x. Then ttIS) is a semialgebraic set 
in R". 


We start by presenting two examples of semialgebraic sets that are relevant 
for our purposes. 

Lemma 1. The set 6„ of stable n x n real matrices (i.e., those with spectral 
radius less than one), when viewed as a subset o/R" , is semialgebraic. 

Proof. An equivalent characterization of stable matrices is via quadratic Lya¬ 
punov functions: 

6„ = {A: 3P = P^ s.t. 

x^Px > 0,yx 0, (3) 

x'^{P — PA)x > 0, Va: yf: 0}. 

Let Tn be the set of n x n matrices A and P that satisfy these constraints; i.e., 
P and P — A^PA both positive definite. The set Tn is semialgebraic (in fact 
basic semialgebraic). One way to see this is to note that the variables a: in ([3]) 
can be eliminated: Since a matrix is positive definite if and only if its n leading 
principal minors are positive, we can describe the set T„ by 2n strict polynomial 
inequalities (in variables that are the entries of P and A). The set ©„ is then 
the projection of T„ onto the space of the A variables. Hence, by Theorem [21 
6„ is semialgebraic. □ 

Lemma 2. The set Vn,d of nonnegative polynomials in n variables and (even) 
degree d is semialgebraic. 

Proof. This is a standard fact in algebra; see e.g. [9]. A polynomial p is by 
definition nonnegative if 


Va:iVa :2 .. .Vx„ p{xi,... ,Xn) > 0. 
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One can apply quantifier elimination techniques (see e.g. m) to eliminate the 
quantified variables (xi,... ,Xn) and obtain a description of Vn,d as a semial- 
gebraic set in terms of the coefficients of p only. Another approach is to resort 
to the representation of nonnegative polynomials as sums of squares of rational 
functions (a la Hilbert’s 17th problem [35]) and note the equivalent description: 

= {p ■ G R[x] of degree d = d{n, d) s.t. 
q is SOS, 
pq is sos}. 

By definition of the sum of squares decomposition, the last two conditions can 
be written as polynomial equations in the coefficients of the polynomials p, q, 
and the coefficients of the new polynomials that are being squared. Hence, the 
result is a semialgebraic set. The set 'Pn,d is the projection of this set onto the 
space of coefficients of p and is therefore semialgebraic. □ 

Unlike the case of stable matrices (Lemma |T|), when we move to switched 
systems defined by even only two matrices, the set of stable systems no longer 
defines a semialgebraic set. This is a result of [33]. The result is stated in 
terms of the joint spectral radius (see |21] for a monograph on the topic) which 
captures the stability of a linear switched systeir@. 

Definition 2. ( \3T^ ) If II • II is any matrix norm, for any integer k, consider 

Pfc(E) := sup{||Ai...Afc||i/'=}. 

AiGE 


r/ie joint spectral radius (JSR) of Ai is 

p(E) = lim pfe(S). 

k—¥oo 


(4) 


The joint spectral radius does not depend on the matrix norm chosen thanks 
to the equivalence between matrix norms in finite dimensional spaces. It is well 
known that the switched system in m is absolutely asymptotically stable if and 
only if p(E) < 1. 

Theorem 3 ( [l^; see also [l^). The 

set of 2x2 matrices Ai, A 2 with p(Ai, A 2 ) < 1 is not semialgebraic. 


Proof. The proof of Kozyakin is established by showing that the family of ma¬ 
trices 

( 1-0 

^ (l-37rt3/2) 


y/l — —t 

0 0 


A2 


( 1-0 


vr^ 

t 


—t 


^ See | 43| for a recent software toolbox for computation/approximation of the Joint Spectral 
Radius. 
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have JSR less than one for t = sin (for k G N large enough), and JSR 
more than one for t = sin (for fc C N large enough). Hence, the stability 
set has an infinite number of disconnected components and therefore cannot be 
semialgebraic. □ 

We now show that by contrast, for any integer d, the set of matrices {Ai, ..., A^} 
that admit a common polynomial Lyapunov function of degree < d is in fact 
semialgebraic. This establishes the result related to the first construction in 
Theorem [T] 

Theorem 4. For any positive integer d, the set of matrices {Ai,...,Am} 
(viewed as a subset o/K™" ) that admit either (i) a polynomial Lyapunov func¬ 
tion of degree < d, or (ii) a polytopic Lyapunov function with < d facets, or 
(Hi) a piecewise quadratic Lyapunov function (inform of max-of-quadratics or 
min-of-quadratics) with < d pieces is semialgebraic. 

Proof. We prove the claim only for polynomial Lyapunov functions. The proof 
of the other claims are very similar. The goal is to show that the set 

'In.m := {Ai,...,Am: 3p := p{x), degree{p) < d, s.t. 

p{x) > 0,\/x 0, . . 

p{x) — p{Aix) > 0,yx 0, 

\/i = 1,..., m} 

is semialgebraic. Our argument is a simple generalization of that of Lemma [T] 
and relies on Lemma [H It follows from the latter lemma that for each d G 
{2,..., d}, the set of polynomials p of degree d and matrices {Ai,..., A^} that 
together satisfy the constraints in ([5]) is semialgebraic. (Note that these are sets 
in the coefficients of p and the entries of the matrices A^, as the variables x are 
being eliminated.) The union of these sets is of course also semialgebraic. By 
the Tarski-Seidenberg theorem (Theorem[2]), once we project the union set onto 
the space of variables in Ai, we still obtain a semialgebraic set. Hence, is 
semialgebraic. □ 

3 Complexity of Lyapunov functions and the finiteness 
property of optimal products 

A set of matrices {Ai,..., Am} satisfies the finiteness property if its JSR is 
achieved as the spectral radius of a finite product; i.e., if 

p(Ai,...,A™) = pi/'=(A<,, ...A,,,), 

for some k and some (ufe,..., ui) G {!,..., m}^. The matrix product Ao-j, ... Ag-^ 
that achieves the JSR is called the optimal product and generates the “worst case 
trajectory” of the switched system in ©■ The finiteness conjecture of [55] (see 
also [IH], where the conjecture is attributed to Pyatnitskii) asserts that all sets 
of matrices have the flniteness property. The conjecture was disproved in 2002 
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by [H] with alternative proofs consequently appearing in [T2], [M], and ^D] . 
In particular, the last reference provided the first explicit counterexample only 
recently. It is currently not known whether all sets of matrices with rational 
entries satisfy the finiteness property ( [H]). 

Gurvits has shown [TB] that if the set of matrices admits a polytopic Lya¬ 
punov function, then the finiteness property holds. The result has been gener¬ 
alized by Lagarias and Wang [2^ to Lyapunov functions that take the form of 
various other norms, including ellipsoidal norms. For example, their result has 
the following implication for quadratic Lyapunov functions. 

Theorem 5 ( [26|1. The finiteness property holds for any set of nx n matrices 
{yli,. . . , of JSR equal to one that share an ellipsoidal norm, i.e., satisfy 
AfPAi :< P for some symmetric positive definite matrix P. Moreover, the 
length of the optimal product is upper bounded by a quantity that depends on n 
and m only. 

In this section, we further generalize this result to “piecewise algebraic” Lya¬ 
punov functions (see definition below). This turns out to include all the classes 
of Lyapunov functions treated in this paper. Note that polynomial Lyapunov 
functions of degree > 4 do not in general define a norm as their sublevel sets 
may very well be non-convex. This is the reason for the need to work with a 
more general version of the result of Lagarias and Wang. 

A minor technical remark before we get to our generalization: In Theo¬ 
rem [5] (and in Theorem |6] below), the assumption is that the JSR is equal to 
one (i.e., the system is neither exponentially stable nor exponentially unsta¬ 
ble). However, the result trivially applies to any positive value of the JSR, by 
homogeneity of this quantity. The condition would instead become that there 
exists a quadratic function V (resp. a piecewise algebraic function V) such that 
Va; G M", VA G {Ai ,..., A^}, V{Ax) < pV{x). 

By a piecewise algebraic Lyapunov function V (of degree d), following the 
nomenclature of [26j . we mean a positive definite, homogeneous, and continu¬ 
ous function (of degree d) whose unit sphere is included in the zero level set 
of a polynomial p{x), with the property that p(0) ^ 0. We refer to p as the 
“associated polynomial” of V. 

We remark that the class of piecewise algebraic functions contains all types of 
Lyapunov functions presented in our introduction section. For example, if V is 
a pointwise minimum of polynomials qi,, qm, then the associated polynomial 
would be p{x) = Y\T=i^t{x) — 1). (Note that the unit sphere of V is the 
union of the unit spheres of the polynomials q^.) Similarly, if F is a pointwise 
maximum of polynomials qi,... ,qm, then the associated polynomial would be 
P{^) = — 1)^. (Note that in this case the unit sphere of V is the 

intersection of the unit spheres of the polynomials qi.) 

We now present the generalized result of interest. The idea of the proof 
follows that of [26] . However, we present here a simplified version, neglecting 
the quantitative estimation of the upper bound. 
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Theorem 6 . Let {Ai ,..., Am} he a set of n x n matrices of JSR equal to one. 
If there exists a piecewise algebraic Lyapunov function V (x) that satisfies 

Vf = 1,..., m, Vx G K”, V{x) — V(Aix) > 0, ( 6 ) 

then {Ai ,..., Am} satisfies the finiteness property. Moreover, the length of the 
optimal product is upper hounded by a quantity that depends on n, m, and d 
only, where d is the degree of the polynomial associated with V{x). 

Proof. Step 1. We start by proving that if p = 1 and V satisfies m, there 
exists a point x* G K", with V{x*) = 1, together with a sequence of matrices 
Ai^,Ai.^,..., such that for all t G N, V(Ai^ .. .Ai^x*) = 1. Clearly, by dH), for 
any x and for any ii,i 2 , V{Ai^ .. < V{x). Now, for any finite t, 

there exists an xt, V(xt) = 1, and a product Ai,. ^ ... ^ such that 

= 1 - 

(Indeed, the contrary would imply that Vx, Vfi ,.. .it, V{Ai^^ ... Ai^x) < 7 < 1, 
which would contradict p = 1.) In fact, by dH), the product has the stronger 
property that Ai^ j/ • ■ • ^ 7,1 = I Vt' < t. 

Finally, by compactness of the set {x : V{x) = I}, and compactness of the set of 
words ii.. .it : t G N under the standard metriqfl, the sequence {xt, (*i, ■ ■ • it)) 
has a subsequence which converges to some (x*, (li, I 2 , ■ • •)) F(x*) = 1. By 

continuity of V{x), we have V{Ai^ ... Ai^x*) = 1 for alH > 1. 

The rest of the proof follows the reasoning of [2^. 0 

Step 2. Now, from the sequence /i, I 2 ■ • ■ obtained above we extract a new 
one ji, j 2 • • • = lAs+i ■ ■ ■ for some s with the supplementary property that there 
are infinitely many 4 G N such that jij2 ■ ■ ■ jt^ = ftk+L^ ■ ■ ■ 3t,,+Lk+tk-i for 
some Lk > 0 '. that is, infinitely many prefices of the sequence turn out to 
reappear later in the sequence. There always exists such a sequence ji,j 2 , ■ ■ ■ 
(see [26l Lemma 3.2]). 

Step 3. Let us denote 

St = {x gM." : V{x) < I, p{Aj^x) =0, ..., p{Ap ... Aj^x) = 0}. (7) 

These sets are nonempty (as they contain x*) and bounded away from zero (as 
p(0) 0). By a standard algebraic geometry argument, the sets St can only 

take finitely many different values (see [551 Lemma 3.3]). They also are such 
that 

yt,St+iCSt. ( 8 ) 

(All the inclusions are non-strict here.) Thus, there exists some tk (as defined 
in step 2 above) such that St,, = S2tk+Lk-i- Plugging this last relation in dZl), 
m, we obtain 

fc +Z;fe — 1 • ■ • -^1 fc+Lfc — 1 ^ ^2tk-\-Lk — l * 

® Take the following distance between two (finite or infinite) words: d{wi,W 2 ) = 1/(2^), 
where I is the index of the first character for which wi{l) ^ W2{1)- The set of (finite and 
infinite) words is compact under this metric. 

^ Our version is however simplified for our purposes and does not effectively compute the 
upper bound in the statement of the theorem. 




3 Complexity of Lyapunov functions and the finiteness property of optimal products 12 


Final Step. Thus, we have a compact set S, bounded away from zero, 
together with a product A = At^^+Lk-i ■ ■ ■ Ai such that AS C S'. As a conse¬ 
quence, for any integer T, (A)^S C S. Since S is compact and bounded away 
from zero, p{A) = 1, and the proof is done. 

It is possible to effectively bound the length tk + Lk of the optimal product 
in the worst case; see [55] for details. □ 

We remark that the theorem above positively answers a previous conjecture 
of ours (see [3] after Theorem 6). We will need a quick lemma: 

Lemma 3. The matrices in 

M. = {Ai,A2} = "I 

satisfy: 

Vx yf 0, |{Ax/||Aa:|| : A G M*}\ = oo, (9) 

where Ai* is the set of all products of matrices in Ai. 

Proof. One can verify that the product A = A 1 A 2 has two real eigenvectors 
with eigenvalues of different modulus. Thus, if x is not an eigenvector of A it 
is clear that |{A*a;/||A‘a;|| : t G N}| = 00 . Now, if X is an eigenvector of A, 
one can check that Aix is not, and we can reiterate the argument on the set 
{A*Aix/\\A*Aix\\:tGN}. □ 

We can now state the results that we were after. 


1 1 
0 1 


1 0 
1 1 


Corollary 1 . Let Ma be the set containing the following two matrices 


with 


Ai 


1 

0 


Ip, A2 — > 


0 

1 


Ip, 


fl 

1' 


1 

o' 

0 

_ 1 

1 

, Ot* 

1 

1 


a* ~ 0.749326546330367557943961948091344672091. .0. 


( 10 ) 


For any positive integer d, there exists a positive number a < a^, such that the 
set of matrices Aia is asymptotically stable under arbitrary switching but does 
not admit a polynomial Lyapunov function of degree < d. 


Proof. Suppose by contradiction that each of these values of a admits a degree d 
or smaller positive homogeneous polynomial Lyapunov function . We suppose 
on top of that that the maximal coefficient of these polynomials 14 , is one 
(scaling the coefficients if needed). Thus, the functions 14, tend to a limit V = 
limc_>a» 14 ,, which is again a homogeneous polynomial function with coefficients 
smaller or equal to one and degree smaller or equal to d because the set of such 

® See 1201 for an expression for the exact value of a*. It is known that violates the 

finiteness property. 
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functions is a compact set. Moreover, it is easy to see that this limit function 
is a valid Lyapunov function for every a < a*. Thus, M.a, must satisfy the 
nonstrict inequalities: 


'iA G M-a ,, V[Ax) < V{x). 


( 11 ) 


Before we conclude, it remains to show that 1C is a valid piecewise algebraic 
Lyapunov function as defined above. The only nontrivial part is that V is 
positive definite, i.e., that there is no a: G R", ||a:|| = 1, such that V{x) = 0. In 
order to prove that, we make use of Lemma |31 The lemma implies that if V has a 
certain x, ||a:|| = 1, such that l^(a:) = 0, by Equation (llip . V has infinitely many 
zeros on the Euclidean sphere, and thus it cannot be a homogeneous polynomial 
of degree at most d with maximum coefficient equal to one. 

We can thus apply Theorem| 6 l which implies that Aia, satisfies the finiteness 
property, a contradiction. □ 

The next corollary is very similar but the matrix family that it presents is 
completely explicit. 

Corollary 2. Consider the matrix family 
(1 - \){Ai,A 2 }, with 


Ai=a'^ 


0 0 
1 0 


, A2 = a 


cos 


2k 


— Sin 


2k 


Sin 

cos 


where 

l<a< (cos^)”^ 

For any positive integer d, there exists a positive integer k such that the set of 
matrices (1 — ^){Ai, A 2 } is asymptotically stable under arbitrary switching but 
does not admit a polynomial Lyapunov function of degree < d. 


Proof. Consider the matrix family {Ai{k),A 2 [k)}. It is shown in that the 
JSR of this family is 1 for all k. On the other hand, the spectral radii of all 
products of length < k are less than one, whereas there is a product of length 
k + 1 which achieves the JSR; i.e., has spectral radius one. As a result, by 
increasing fc, the length of the optimal product blows up. Hence, by Theorem| 6 l 
the degree of a polynomial Lyapunov function (applied to the asymptotically 
stable family (I — ■i){Ai(fc), A 2 [k)}) must blow up as well. □ 

In Section O we test the result of the above corollary by searching for poly¬ 
nomial Lyapunov functions using semidefinite programming. 


4 Complexity of Lyapunov functions and undecidability 

In this section, we show that our statements on lack of upper bounds on com¬ 
plexity of Lyapunov functions also follow in a straightforward manner from 
undecidability results. Compared to the results of the previous sections, the 
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new statements are weaker in some sense and stronger in some other. They are 
weaker in that the statements are non-constructive. However, they imply the 
stronger statement that the complexity of Lyapunov functions (e.g., degree or 
number of pieces) cannot be upper bounded, not only as a function of n and m, 
but also as a computable function of n, m, and the entries of the matrices in E 
(Corollary 0]). In addition to this, we can further establish that the same state¬ 
ments are true for very simple and restricted classes of matrices whose entries 
take two different values only (see Theorem ITUl) . 

Theorem 7. For any positive integer d, there are families of matrices of size 
47 X 47 that are asymptotically stable under arbitrary switching but do not admit 
(i) a polynomial Lyapunov function of degree d, or (ii) a polytopic Lyapunov 
function with d facets, or (Hi) a max-of-quadratics Lyapunov function with d 
pieces, or (iv) a min-of-quadratics Lyapunov function with d pieces. 

The main ingredient in the proof is the following undecidability theorem, 
which is stated in terms of the JSR of a set of matrices. 

Theorem 8. ( UlVil.Sf ) The problem of determining, given a set of matrices E, 
if p(E) < 1 fs Turing-undecidable. This result remains true even if E contains 
only two matrices with nonnegative rational entries of size 47 x 47. 

We now show that this result implies Theorem [T] The main ingredient 
is Tarski’s quantifier elimination theory, which gives a finite time procedure for 
checking certain quantified polynomial inequalities. The rest is a technical trans¬ 
formation of the problem “p < 1 ?” to the existence of a degree d polynomial 
Lyapunov function. 

Proof of theorem We suppose by contradiction that every AAS set of matrices 
of size 47 x 47 has a polynomial Lyapunov function of degree at most d, for some 
even natural number d. We claim that this implies the algorithmic decidability 
of the question “p < I?”. 

Indeed, by homogeneity of the JSR, we have 

p(E)<IoVee(0,l), p((I-e)E)<l. 

Now, by the hypothesis above, the last statement should be equivalent to the 
existence of a polynomial Lyapunov function of degree less than d for (1 — e)E, 
and thus we can rephrase the question p(E) < 1 as follows. (In what follows, 
Vnci of polynomials in n variables and (even) degree d that are positive 

definite. This set is semialgebraic. See Lemma [2] for a similar statement.) 

Ve S (0,1), 3p(-) e 7^+^ such that 
\/x G \ {0}, VA G E, p((I — e)Ax) < p{x). 

The above assertion is algorithmically decidable via Tarski’s quantifier elim¬ 
ination theory [41]. Thus, it allows to decide whether p(E) < 1, contradicting 
Theorem |5| 
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The proof for the polytopic Lyapunov function and max/min-of-quadratics 
goes exactly the same way, by noticing that once their number of components 
is fixed, one can decide their existence by quantifier elimination as well. □ 

In [13], the authors note that Theorem |S| implies the following result: 

Corollary 3. ( fTS^) There is no effectively computable functioi^ which 

takes an arbitrary set of matrices with rational entries E, and returns in finite 
time a natural number such that 

pCEi) = max maxp(A). 


The same corollary can be derived concerning the degree of a Lyapunov 
function. 

Corollary 4. There is no effectively computable function (i(E), which takes 
an arbitrary set of matrices with rational entries E, and returns in finite time 
a natural number such that if p(E) < 1, there exists a polynomial Lyapunov 
function of degree less than d. 

Next, we show a similar result, which does not focus on the fixed size of 
the matrices in the family, but somehow on the complexity of the real numbers 
defining the entries of the matrices. Namely, we show that such negative results 
also hold essentially for sets of binary matrices (that is, matrices with only 0/1 
entries). In fact, the very question p < 1 is easy to answer in this case (see [23]L 
so, one cannot hope to have strong negative results stated in terms of binary 
matrices. However, it turns out that for an arbitrary integer K the question 
p < K for binary matrices is as hard as the question p < 1 for rational matrices. 
More precisely, we have the following theorem: 

Theorem 9. ( i2Sf ) Given a set of m nonnegative rational matrices E, it is 
possible to build a set of m binary matrices E' (possibly of larger dimension), 
together with a natural number K such that for any product A = S E*, 

the corresponding product ... A'^ G E'* has numerical values in each of its 
entries that are either exactly egual to zero, or to some entry in the product A 
multiplied by K*. Moreover, for any entry in the product A, there is an entry in 
the product A' with the same value multiplied by K*. 

Thus, the above theorem provides a way to encode any set of rational ma¬ 
trices in a set of binary matrices (in the sense that (i) the joint spectral radius 
is the same, up to a known factor K, and (ii) that there is a perfect correspon¬ 
dence between products of matrices in E and in E'). Theorem [8| together with 
Theorem ini allows us to prove another negative result on the degree of Lyapunov 
functions restricted to matrices with entries all equal to a same number 1/K, 
K G Q. Remark that the fact that the parameter K G Q, has unbounded denom¬ 
inator and numerator is unavoidable in such an undecidability theorem, since 


® See ( [13]) for a definition. 
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for bounded values, there is a finite number of matrices with all entries in the 
range, and this would rule out a result as the one in the theorem below. 

Theorem 10. There is no computable function d : N x Q —> N such that for any 
set of matrices of dimension n with entries all taking values in the set {0,itr} 
for some if S Q, the set is A AS if and only if there exists a (strict) polynomial 
Lyapunov function of degree d{n, K). 

Proof. Given a set of rational matrices of size 47 x 47, Theorem 0 allows us 
to build a set of binary n' x n' matrices S' such that p(S') = Kp{Yi). Thus, 
p(Yj'!K) = piYAj, and the existence of strict polynomial Lyapunov functions for 
AAS n' X n' matrices would again imply decidability of the question p < 1 for 
Yf jK., which again contradicts Theorem [ 8 ] □ 

Roughly speaking, our last theorem shows that in higher dimensions, “bad” 
families of matrices that necessitate arbitrarily complex Lyapunov functions can 
have very simple and structured entries. 

5 Numerical validation 

In this section, we test the result of Theorem [T] through a numerical experiment. 
This will also give us a feel for how quickly the complexity parameter of the 
Lyapunov function can blow up in practice. 

As a representative example, we choose the matrices of Lagarias and Wang 
from |26] . parameterized by a single integer k. For each k > 0, we know the 
matrices have JSR< 1, but we try to find the smallest degree sum of squares 
Lyapunov function that can prove this. A nice feature of this example is that 
since the matrices are 2x2, the Lyapunov functions and their increments will 
be bivariate forms. It is well-known that all nonnegative bivariate forms are 
sums of squares (see, e.g., [36]). Hence, the minimum degree of an sos Lay- 
punov function is also exactly the minimum degree of a polynomial Lyapunov 
function. However, working with the sos formulation allows us to to compute 
this minimum degree using semidefinite programming. See, e.g., [32], [ 6 ] for 
an explanation of how this is done. Our SDP solver throughout the example 
is SeDuMi |40j . It would be worthwhile to perform similar experiments with 
min/max-of-quadratics or polytopic Lyapunov functions, etc. 

The input matrices to our problem are the following 

cos t sin 

-Sint cos^J’ 

where 

a = ^(1 + (cos^)"^), 

and fc is a parameter. Our results are as follows: 


Ai — (1 — 


IQSfc' 


0 0 
1 0 


A 2 = (1 — 


IQSfc' 
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k 

2 

3 

4 

5 

6 

7 

d 

6 

10 

14 

16 

20 

24 


This table is stating, for A: = 7 for example, that the matrices are stable 
under arbitrary switching and yet this fact cannot be proven with a polynomial 
Lyapunov function of degree less than 24. This trend keeps growing as expected 
until we run in numerical issues with accuracy of the SDP solver. We believe 
this example can now serve as a parametric benchmark problem for comparing 
the power of optimization-based techniques for proving stability of switched 
systems. 

6 Conclusion 

In this paper, we leveraged results related to non-algebraicity, undecidability, 
and the finiteness property of the joint spectral radius to demonstrate that com¬ 
monly used Lyapunov functions for switched linear systems can be arbitrarily 
complex, even in fixed dimension, or for matrices with lots of structure. 

If these negative results are bad news for the practitioner, it is worth men¬ 
tioning that in practice the different Lyapunov functions often have comple¬ 
mentary performance. So while there certainly exist instances which make all 
methods fail (as we have shown), one can hope that in practice, at least one of 
the different Lyapunov methods would be able to certify stability. In light of 
this, we believe it is important to (i) understand systematically how the different 
methods compare to each other, and (ii) identify subclasses of matrices that if 
stable, are guaranteed to admit “simple” Lyapunov functions. While the latter 
research objective has been reasonably achieved for quadratic Lyapunov func¬ 
tions, results of similar nature are lacking for even slightly more complicated 
Lyapunov functions (say, polynomials of degree 4, or piecewise quadratics with 
2 pieces). 
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